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ABSTRACT 

We present a new, approximate method for modelling the acceleration and collimation 
of relativistic jets in the presence of gravity. This method is self-similar throughout the 
computational domain where gravitational effects are negligible and, where significant, 
self-similar within a flux tube. These solutions are applicable to jets launched from 
a small region (e.g., near the inner edge of an accretion disk). As implied by earlier 
work, the flow can converge onto the rotation axis, potentially creating a collimation 
shock. 

In this first version of the method, we derive the gravitational contribution to the 
relativistic equations by analogy with non-relativistic flow. 

This approach captures the relativistic kinetic gravitational mass of the flowing 
plasma, but not that due to internal thermal and magnetic energies. A more sophis- 
ticated treatment, derived from the basic general relativistic magnctohydrodynamical 
equations, is currently being developed. 

Here we present an initial exploration of parameter space, describing the effects the 
model parameters have on flow solutions and the location of the collimation shock. 
These results provide the groundwork for new, semi-analytic models of relativistic 
jets which can constrain conditions near the black hole by fitting the jet break seen 
increasingly in X-ray binaries. 

Key words: acceleration of particles - ISM: jets and outflows - methods: analytical 
- MHD. 



1 INTRODUCTION 

Despite decades of study, astro-physical jets are still an enig- 
matic component of our universe. Understanding the details 
of their physics would advance fields as diverse as star for- 
mation, gamma-ray bursts (GRBs), galaxy evolution and 
cluster dynamics. Jets in all of these objects probably share 
three basic ingredients: a source of energy and particles such 
as an accretion flow, magnetic fields created or carried in- 
ward by the flow, and rotation. However, the exact details 
of how these elements combine to launch and collimate jets 
in the observed systems is not at all clear. 

Black holes are ideal test cases to understand jet 
physics because of their enormous range in mass, allow- 
ing us to study jet dynamics over an equally large range 
of time-scales. Accreting stellar mass black holes in a bi- 
nary with a companion star (BHBs) have jets that can un- 
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dergo an entire launch and quench cycle on time-scales of 
months, which is related at least in part to the accretion 
state of the disc jEsin et al.lll997l; lMeierll200ll ; iFender et all 
12004 iMcClintock fc RemiHardll2006t ). In their supermassive 



cousins, active galactic nuclei (AGN), such transitions would 
be expected to occur on time-scales of millions to a billion 
times longer. To assess the extent to which such a mass- 
scaling operates, we must have a better understanding of 
which elements of jet physics persist despite a dramatic 
range of environment and scales. 

One key feature of jets regardless of source seems to 
be their ability to accelerate particles to highly relativistic 
energy. We observe the outcome of these accelerated par- 
ticles direct ly via optically thin sy nchrotron emis sion in, 
e.g. AGN (|Marscher fc Gearl Il985l ) and in BHBs (|Fenderl 
l200lf ). When the jets are compact enough to become self- 
absorbed, the superposition of many synchrotron-emitting 
components along the length of the jet leads to a flat or 
slightly inverted total spectrum l|Blandford fc Konigllll979h 



© 2012 RAS 



2 Peter Polko, David L. Meier and Sera Markoff 



until the point where the jet emission becomes entirely op- 
tically thin. At this point the spectral index will steepen 
dramatically to a ~ 0.5 — 0.8 typically, where a is defined 
such that F v oc v~ a . This break in the spectrum can be 
physically associated with the most compact region in the 
jets where particle acceleration is present. In low- lumino sity 
AGN jets this break typically occurs in the GHz range (|Hol 
1 19991 ). and if jets are self-similar across the black hole mass 
range, in BHBs such a break would be predicted to oc- 
cur in the infrared (IR ) bands (jMarkoff et all l200ll . 120031 ; 
iHeinz fc SunvaevlbOOsh . Because the spectral break also in- 
dicates the total radiative power of the jets, and plays a 
key role in driving mass-scaling predicti ons such as the Fun- 
damental Plane of black hole accretio n (|Merloni et al.ll2003l ; 
iFalcke et al]|2004l : |p"lotkin et al.ll2012h . it is a pivotal param- 
eter linking the inner jets (and assumedly conditions near 
their launch point) to the outer jets associated with strong 
radio through IR emission. Deriving this link and overall 
jets structure in an MHD-consistent way is the focus of this 
paper. 

Constraints on the location of the jet break can come 
from both observations as well as from spectral fitting. Be- 
cause the optical/IR (OIR) bands are often dominated by a 
stellar companion or the accretion disc, the optically thick- 
to-thin break ha s been observed explicitly in only one BHB 
so far, G X 339-4 ICorbel fc Fenderl2002l ; lG"andhi et alj|201ll ; 
see also IMigliari et alj l200rj for a similar break for neu- 
tron stars), and has been indirectly co nstrained for Cyg X - 
1 based on recent MIR observations jRahoui et al.l [201lf ). 
These observations confirm the theoretically predicted break 
location, and the recent boom in multiwavelength monitor- 
ing of many BHBs is leading to new studies of the break and 
its scaling with luminosity (Russell et al., in prep.). The lo- 
cation in frequency of the break can also be rather tightly 
constrained by fitting multiple spectra of BHBs and AGN 
in states associated with compact jets, which has the ad- 
vantage of also associating a size scale for the region, and 
distance from the black hole. A series of works have con- 
sistently found the region of the jet where particle accel- 
eration initiates seems to be offset from the black hole at 
distances of ~ 10-1000 r g , depending on sou rce luminosity, 
wher e r E is t he gr a vitational radius, GM/c 2 (iMarkoff et al. 



2001 



120031 . 

iMaii 



2005, 



i. l2008l ; IMigliari et all 120071 ; iGallo et al 
20071 ; iMaitra et alj|2009i ) . Accordingly, this zone also should 
be associated with an offset in the start of the synchrotron 
core from the black hole in direct imaging. Such a measure- 
ment is very difficult to do because of the spatial resolution 
required. At lea st for the jet in M8 7, a relatively nearby AGN 
(17.0±0.3 Mpc. lTonrv et alj|2001) with a large supermassive 



black hole (6.4 ± 0.5 x 10 9 Mm . iGebhardt fc Thomasll2009l . 



although note this value is twice the m ass found in previ- 
ous studies), the offset is ~ 30-100 r g |junor et al.l 1 19991 : 
I Walker et al]|2008l ; Irlada et al.l l201lh . Interestingly, recent 
astrometry for M87 has also shown that the location of the 
radio core with res pect to the black ho le is extremely stable 
over several years (|Asada et alj|201ll ). A natural explana- 
tion for such a region for the start of particle acceleration 
would be a shock where continual diffusive acceleration oc- 
curs (e.g. Hil Il978l ; iDrurvl Il983l ) . Because particles need 
to be continually accelerated, the shock should be a steady 
feature within the compact jets. 

In self-similar, axisymmetric MHD flows, the bulk flows 



can be accelerated through several singular points in order: 
the modified slow point (MSP), t he Alfven point (AP), and 
the modified fast point (MFP) (|Blandford fc Payne! jl982j 
hereafter BP82). At the MFP, the velocity of the flow greatly 
exceeds the fast magnetosonic speed, which is the fastest 
that signals in the jet can travel. Thus beyond the MFP the 
flow is causally disconnected from the flow closer to the cen- 
tral object, allowing disruptions such as shocks to form. A 
stable location for the MFP in a given flow would therefore 
provide a natural explanation for a steady feature associated 
with particle acceleration. This location could, however, be 
a function of mass and accretion power and vary from source 
to source within a certain range. In the models mentioned 
above, the shock location is a fitted free parameter that sug- 
gestively falls into the same range for a variety of sources. 
However, to understand if this is a physically meaningful 
result, we should ideally be able to show that this location 
corresponds to a feature that can be derived self-consistently 
within a theoretical framework, and tied to conditions at the 
jet launch point. The earlier mo des are based on a hydro- 
dynamical (HP) veloci ty profile (|Falcke fc Biermannlll995l ; 
IFalcke fc Markofj|2000h , and thus do not invoke the role of 
magnetic fields other than as a global parameter. However, 
if we move to a MHD framework, we can derive the exact 
location of the MFP (and assumedly the shock) based on 
boundary conditions defined by the launch point of the jets, 
and then test if the theory provides the correct scalings as 
observed, and as seem to be driving the Fundamental Plane. 

In this paper we take a semi-analytical approach to de- 
rive new flow solutions (a new jet model) with the explicit 
aim of exploring the conditions for formation of the MFP, by 
solving the equations of relativistic MHD under the assump- 
tions of self-similarity and axisymmetry. These assumptions 
greatly simplify the equations describing the jets, and are 
consistent with the results of numerical MHD simulations, 
which almost always produce a magnetic field line geome- 
try near the launch point that is remarkab ly self-similar an d 
axisymmetric (see, e.g., fig. 2 and fig. 11 in lMcKinney||20o"r3 ) . 

We need to make one further approximation to con- 
struct our solutions: in order to link the properties of the 
MFP to the conditions at the base of the jet, we need to 
solve for the AP and MSP in a regime which could poten- 
tially be very close to the black hole where gravity becomes 
important. The MSP is a natural point to associate with the 
launch region of the jets, similar to the sonic nozzle in hydro- 
dynamical flows, and is thus also the place where conditions 
could eventually be matched to the accretion inflow, such as 



(RIAF; iNaravan & Yil Il994 


; iBlandford & Begelmanl Il999l: 


iQuataert & Gruzinovll2000l: 


Merloni & Fabian 120021). How- 



ever, gravity will play a significant role so close to the black 
hole; therefore, it also needs to be included in order to 
correctly derive and cross the MSP. Until now no semi- 
analytical formalism has been developed to describe a rela- 
tivistic flow passing through all three singular points, because 
in a relativistic framework gravity is not self-similar. There- 
fore in this paper, in order to accomplish this connection, we 
apply a bridging solution between our p revious relativistic 
self-similar flow model without gravity (|Polko et all |2010| . 
hereafter PMM 10, based on the framewo rk for self-similar 
flow laid out in IVlahakis fc Konigll 120031 , hereafter VK03), 
which is also valid in the non-relativistic limit (again without 
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gravity), and a non-relativistic flow model that is valid close 
to the black hole and that does include the effects of gravity 
self-consistently l|Vlahakis et al.l l200d . hereafter VTSTOO). 
The transition from the non-relativistic flow with gravity 
to our relativistic flow without gravity from PMM10 occurs 
sufficiently far from the black hole that gravitational effects 
are no longer important. The bridging of these two solu- 
tions will not describe all possible gravitational fields and 
velocities, but the range of solutions is wide enough to be 
applicable to many observed sources, which is the goal of this 
project. Specifically, within a cylindrical region constituting 
a "flux tube", that can be associated with the observable 
"sheath" of compact jets, we demonstrate that the devia- 
tions from self-similarity when gravity is present are very 
small. Within this scenario, we present relativistic MHD so- 
lutions that cross all three singular points, allowing us to 
derive the location of the shock self-consistently from bound- 
ary conditions at the base of the jet and to begin to test how 
well such a model fares against physical sources. 

Gravity is treated in the weak (Newtonian) limit, and 
in this paper we have opted for a Paczyhsky-Wiita (pseudo- 
Newtonian) augmentation since some of our solutions ap- 
proach the Schwarzschild radius, where general relativistic 
effects begin to be important. No other GR correction terms, 
nor black hole spin, is treated at this time. 

In §2 we present details of the bridging formalism. In §3 
we use the new formalism to obtain solutions crossing the 
MSP, AP, and MFP and describe their properties. We also 
show the effects that changing the model parameters have 
on the location of the MSP and MFP. In §4 we present an 
in-depth discussion of our model and results, and we draw 
our conclusions in §5. 



2 METHOD 

In order to tie the conditions at the start of the particle 
acceleration region to the conditions at the base of the jet, 
we need to have a full description of a hot, relativistic flow. 
Because close to the central object gravity is as important 
as temperature and the magnetic field strength, the gravi- 
tational potential must be taken into account in order to be 
able to describe the jet in this region. Because this region is 
generally subrelativistic or mildly relativistic, we start with 
a solution that includes gravity in a non-relativistic flow. Far 
from the black hole the flow will be relativistic, but grav- 
ity will no longer be important in the solution, so we can 
use the equations derived in PMM10 to describe this flow. 
We then seek solutions that satisfy both the relativistic flow 
without gravity and the non-relativistic flow with gravity at 
opposite ends of the jet. 

2.1 A physical description of the flow 

In the framework of self-similar relativistic MHD, we find 
that jets in our solutions are accelerated in the following 
way: starting from the base of the flow, the initial bulk ac- 
celeration is provided through thermal pressure gradients by 
means of sound waves by the very hot particles surround- 
ing the central object. When the flow exceeds the sound 
speed, this mechanism becomes inefficient, and acceleration 
due to the centrifugal force of the rotating magnetic field by 



means of Alfven waves takes over. After the flow surpasses 
the Alfven speed, magnetic pressure becomes the dominant 
mechanism. When the flow velocity increases beyond the fast 
magnetosonic speed, the final boost is given by the pinching 
of the magnetic field. If at this point the flow overcollimates, 
a shock may form, which can accelerate the particles into the 
observed power-law distributions, thus constituting the start 
of the particle acceleration region. 

In prior work (PMM10) we were unable to probe the 
region where the initial bulk acceleration occurs due to the 
absence of gravity in our formalism. However, the singu- 
lar MSP can be produced only by the inward gravitational 
force balancing the outward magnetic and thermal pressure 
forces. Thus by including gravity we can describe the jet 
from very close to the central object to the first instance of 
overcollimation, providing the connection between the con- 
ditions at these two points. We will discuss how we deal with 
the issues regarding self-similarity in S0] 

2.2 A new solution technique: the C°°-continuous 
bridging method 

The matching of two flow solutions valid in different regimes 
of parameter or physical space is a common technique in 
physics. Some techniques simply solve the two different 
equations in the different regimes and then match the so- 
lutions at a single (and arbitrary) point. Such bridging solu- 
tions are only C°-continuous, meaning its derivative is not a 
continuous function. Other methods may involve spline fit- 
ting the two solutions together in a finite bridging region. In 
this case, the two solutions are valid in each of their respec- 
tive regimes, but in the bridging region neither is strictly 
valid. Instead, the spline fit is an approximation of the tran- 
sition between the two solutions, and it is C 1 -continuous or 
better. 

Our method for bridging the non-relativistic VTSTOO 
MHD wind equation with gravity and the relativistic VK03 
one without gravity is similar to the above approaches. That 
is, when gravitational forces are important and the flow non- 
relativistic, we obtain the VTSTOO MHD solution. However, 
when gravity is no longer important, the flow has the charac- 
ter of the VK03 relativistic solutions, whether it is relativis- 
tic or not. The difference here is that, because of the nearly- 
identical structures of the VTSTOO and VK03 equations, 
we can perform this task with the construction of a single 
wind equation. The method, then, will be C^-continuous, 
with all derivatives being continuous functions - a distinct 
advantage when one is dealing with equations which have 
singular points through which the solution must pass. 

2.3 The basic C°°-continuous method 

We now give a brief overview of our previous work in PMM10 
and describe how we here extend it by bridging it with a non- 
relativistic formulation including gravity valid close to the 
black hole. 

In PMM10 we combined the Bernoulli equation (de- 
scribing the forces along a field line) and the transfield equa- 
tion (describing the forces perpendicular to a field line), to 
construct a wind equation, which describes the bulk acceler- 
ation of the flow and fully specifies the jet geometry. Using 



© 2012 RAS, MNRAS 000,rflfl3l 



4 Peter Polko, David L. Meier and Sera Markoff 



this single differential equation, it was possible to obtain jet 
solutions that crossed an MFP at a finite distance from the 
origin. The terms in the wind equation thus obtained using 
the formalism of VK03 turned out to be almost identical 
to the terms in the wind equation given in VTSTOO. Apart 
from relativistic corrections to all the corresponding terms, 
the only difference is the gravity term in the VTSTOO wind 
equation, which is non-relativistic. If we modify the VK03 
wind equation to include gravity, our combined wind equa- 
tion will reduce exactly to the formalism of VTSTOO when 
relativistic effects apart from gravity are negligible as is the 
case close to the black hole, while it also reduces exactly 
to the formalism of VK03 when gravity becomes negligible 
as is the case far away from the black hole. Between these 
two regimes the wind equation describes the bridging regime 
where relativistic effects may begin to dominate the gravi- 
tational effects. Since all these regimes are described by one 
continuous equation, the resulting solution will be a smooth 
transition between the two regimes. We prefer this approach 
over matching the solutions of the two regimes at an arbi- 
trary location in the jet, while trying to satisfy any number 
of continuity conditions. 

Since the notation in VTSTOO differs from VK03, we 
use equations (8)-(12) in VTSTOO and equations (24) in 
VK03 to translate the notation to that used in the latter. 
This conversion also provides the relativistic corrections to 
the gravitational mass. Using this relativistic velocity in the 
definition of the gravity term in VTSTOO yields a relativistic 
form of the gravity term, allowing us to include it into the 
VK03 wind equation. The term by term comparison of the 
two wind equations ensures that we use the proper scaling. 
Since some of our solutions get very close to the central 
black hole, we adapte d the gravity term to inclu de a pseudo- 
Newtonian potential (|Paczvriskv fc Wiitd[l98ch . 

After these steps (see appendix[X]) the gravity term has 
the following form: 



H 2 x\ (1 - M 2 



F 2 



el (1 



M 2 — x 2 )' 2 



c ujf^G 



QM sin(0) 



(1) 



where the subscript A denotes values at the Alfven point and 
/ic 2 is the total energy-to-mass flux ratio, xa is the cylin- 
drical radius distance of the AP scaled to the light cylinder 
radius, F controls the current distribution, <tm is the mag- 
netisation parameter, M is the Alfvenic Mach number, x is 
the cylindrical radius distance scaled to the light cylinder 
radius, c is the speed of light, uja is the cylindrical radius 
distance to the AP, G is the cylindrical radius distance scaled 
to the AP, Q is the gravitational constant, M is the mass of 
the black hole and 6 is the spherical polar angle (see Table 
[2] for an overview of all model parameters). 

The first two fractions consist purely of constants along 
a given field line, providing the overall scaling. The third 
fraction is proportional to (7^) 2 (see I Alf) . where 7 is 
the Lorentz factor and £c 2 is the specific relativistic en- 
thalpy (per baryon mass). The last fraction corresponds to 
(r/r g — 2) _ , with r the spherical radius and r g the gravita- 
tional radius, and it is this term that represents the pseudo- 
Newtonian potential. 

Several equations besides the wind equation are needed 
to fully determine a solution, so we have made sure that we 
take into account all equations in which the gravity term 
appears. As explained above, gravity, as described by equa- 



tion |T}, shows up in the numerator of the wind equation. 
By evaluating the wind equation at the AP, we obtain the 
Alfven regularity condition (ARC, see equation IB1[) . which 
also depends on gravity. The ARC allows us to calculate the 
slope of M 2 at the AP, pa, one of the initial parameters of 
the integration. The gravity term, as it appears on the right 
hand side of the ARC in VK03, is given by: 



l-x 2 A 



C ZUa 



(2) 



GM sin(<9 A ) 

Another location in which the gravity term may appear 
is in the Bernoulli equation evaluated at the AP. This equa- 
tion is used to determine fj,, another initial parameter of the 
integration (see equation IB2[) . However, by evaluating the 
Bernoulli equation of VTSTOO at the AP, the gravity term 
in that equation vanishes there. So equation (B5) in VK03 is 
still valid (see appendix[C]). The transfield equation is incor- 
porated into the wind equation, but not used independently 
in our calculations, so its dependence on gravity is irrelevant 
here. 

By making these changes we have been able to find 
solutions that can pass through all three singular points of 
a hot, relativistic MHD flow, extending all models so far 
published (for a qualitative overview see Table [TJ. 



2.4 The effects of including gravity on the 
solutions 

In every instance of gravity the ratio vja/M appears. Thus 
for a central object that is twice as massive, if the AP is 
moved twice as far from the axis of symmetry, the solution 
remains unchanged. This result is a direct consequence of 
mass-scaling, since all properties are expressed in gravita- 
tional radii. Regulating the gravitational force exerted by 
the compact object can therefore be achieved by changing 
either M or ■cua and, apart from the overall physical size of 
the system, these two actions are equivalent. Increasing the 
effect of gravity can thus be thought of either as increas- 
ing the mass of the central object (increasing the reach of 
the gravitational well), or as decreasing the radius of the 
AP (moving the Alfven point deeper into the gravitational 
well). 

One way to show these effects, would be by starting 
with a singular solution containing both an MFP and an 
MSP, slowly increasing zua to decrease the effect of gravity, 
while keeping solutions with an MFP and MSP. Unfortu- 
nately, these solutions do not smoothly transition into solu- 
tions that only have an MFP. It is therefore impossible to 
directly compare these new solutions crossing all three sin- 
gular points to solutions without an MSP and we are left 
noting the effects within these new solutions. 

The significance of having a solution with an MSP is 
that only then it is possible to connect conditions at the 
MFP to conditions very close to the central object. Even 
though the extension of a field line in linear distance may 
not be much, we can be sure that the results in the sub- 
Alfvenic regime have meaning by satisfying a strict bound- 
ary condition. We equate the MSP with the launch point of 
the jet and fit the conditions at the MSP to the accretion 
flow. Hence, extending the solutions to the MSP allows a 
link between the accretion flow and the start of the particle 
acceleration region, which we will use in later work. 
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Table 1. Overview of self-similar, axisymmetric MHD models, classified whether they include gravity, allow for a warm flow, allow 
relativistic velocities, and cross the MSP, AP and MFP, respectively. 



Paper 


Gravity 


Warm 


Relativistic 


MSP 


AP 


MFP 


Blandford & Pavne (1982) 


/ 








/ 




Li et al. (1992) 






/ 




/ 




Vlahakis et al. (20001 


/ 


/ 




/ 


/ 


/ 


Vlahakis & Koniel (20031 




/ 


/ 




/ 




Polko et al. (20101 




/ 


/ 




/ 


/ 


Polko et al. (this work) 


/ 


/ 


/ 


/ 


/ 


/ 



Table 2. List of model parameters 



Fundamental parameters 


F 


Sets the radial dependence of the magnetic field strength 


r 


Adiabatic index 


x A 


Cylindrical radius in terms of the light cylinder radius 


CM 


Magnetisation parameter 


1 


Dimensionlcss adiabatic coefficient 


TXA 


Cylindrical radius 


M 


Mass of the central object 


Derived parameters 


0A 


Angle with respect to the axis of symmetry 


IpA 


Angle of the field line with respect to the disc 


PA 


Derivative of the Alfvenic Mach number squared w.r.t. 9 


<?A 


Magnetisation func. (Poynting-to-matter energy flux ratio) 


fie 2 


Total energy-to-mass flux ratio 




Specific (per baryon mass) relativistic enthalpy 



Notes: The subscript A means the value of the variable at the AP. 
For a complete description of the parameters, please see PMM10. 



Compared with the Paczyrisky-Wiita potential, the 
Newtonian potential has a weaker gravitational force at the 
same radius. As it is the gravitational force that balances all 
other forces at the MSP, using a Newtonian potential would 
move the MSP closer to the black hole, possibly beyond the 
event horizon for certain solutions. To avoid this, we adopt 
the Paczyhsky-Wiita potential to approximate the gravita- 
tional potential close to the black hole. The adaptation from 
a Newtonian potential is straightforward (see appendix |A| . 

2.5 Model parameters 

By including gravity in the equations, two additional pa- 
rameters need to be specified before an integration. These 
are the mass of the central object, M, and the cylindrical 
radius of the AP, zua- As mentioned above, only the ratio of 
these parameters appears in the equations, so these two pa- 
rameters are linearly dependent and solutions with the same 
ratio are indistinguishable if all other parameters are fixed. 
Since a solution is given in gravitational radii, this means 
that any solution can be scaled up or down to any physical 
size. For our initial solution we have chosen to set M to the 
mass of a typical stellar mass black hole of 10 Mq (although 
this value will eventually be set to the deduced mass of an 
observed black hole system) and to set vja to 5 r s . We will 
look at the effects of mass scaling at a later time. 



The parameters F and Y are set to a single value for 
now, but can be changed if necessary. For example, since F 
influences the geometry of the jet, it also affects the collima- 
tion. By allowing it to vary, we will have the freedom to ad- 
just the degree of collimation to fit the broad-band data on 
our sources. Due to the self-similar nature of the equations, 
the height of the MFP will depend on the initial radius at the 
disc. Because we want a representative location, we choose 
to anchor the field line in the most active part of the accre- 
tion disc close to the central object. Within the innermost 
stable circular orbit (ISCO) the accreting matter can no 
longer form a disc and will move towards the central object 
as a radiatively inefficient flow. Therefore we set F = 0.75, 
which corresponds to a wide range of accretion disc models, 
including the Blandford-P ayne (BP82), Shakura-Sunyaev 
jShakura fc Sunvaevl 1973T I . and radiatively-inefficient accre- 
tion flow (RIAF) models. We also set Y — 5/3, appropriate 
for non-radiation-dominated jets in the hard state, where 
the radiation and the particles do not behave as a single 
fluid. 

The choice of parameters to fit the singular points will 
affect the solution we find. Were we to use 9a and ipA to 
fit for the MFP and MSP positions, the angle of the Alfven 
point and the slope of the field line there would be varied to 
find a singular solution. If, on the other hand, x\ and q were 
used to fit for the MFP and MSP, the light cylinder radius 
(and hence angular velocity) and dimensionless adiabatic 
coefficient would be varied instead. This last combination 
most closely approaches self-similarity out of all the pairs of 
parameters that we have studied. Therefore, in this work we 
will use x\ and q as fitting parameters. It is important to 
note that by choosing different fitting parameters only the 
way we move through parameter space changes, not the set 
of solutions themselves. 



3 RESULTS 

3.1 First solution 

Starting with the warm (initial £ ~ 4.7) and fast (final 
7 ~ 10) solution c in PMM10, we use our new method to 
explore parameter space to determine a solution with an 
MSP. The parameter values of this new solution are very 
close to those of solution c (see Table [3J , establishing the 
relative ease with which solutions can be found. Fig.[T]shows 
the numerator and denominator of the wind equation, and 
their ratio corresponding to the derivative of square of the 
Alfvenic Mach number with respect to the poloidal spherical 
angle 6. From this ratio it is clear we have indeed found a 
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Table 3. Parameters of solutions 
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1 The parameters given are for solution c in PMM10. This is the solution we will compare our new result with. The values for the first 
seven parameters (F through M) of the solution in this work are exact, for the last five {x 2 A through ^t) they are rounded off. Because 
singular solutions require high precision, the rounded-off numbers are given with six significant digits. 



smooth crossing. In comparison with solution c, the AP is 
still located at 8 = 0.87 rad or 50° from the axis of sym- 
metry, but the MFP has moved outward to 9 — 0.029 rad 
or 1.65°, while the MSP (only ostensible in solution c) has 
moved inward to 6 = 0.92 rad or 52.6°. This angle corre- 
sponds to a spherical radius of 5.96 r g , which is just within 
the ISCO. 

The values of the velocity, the Lorentz factor, the mag- 
netic and electric field strength, the density and the pressure 
along a reference field line are given in Fig. [2] The poloidal 
velocity monotonically increases from 0.067 c, through 0.11 c 
at the MSP, to very close to the speed of light at the MFP. 
The toroidal velocity starts off positive, turns negative af- 
ter the AP, and returns to positive before the MFP. The 
Lorentz factor slowly increases from 1.17 at the MSP to a 
final value of 12.3. 

The poloidal and (negative) toroidal magnetic field 
strengths first increase in magnitude and then decrease, both 
peaking around the AP. The electric field has the same be- 
haviour. The density and pressure both drop monotonically 
as is expected for an expanding jet. Just beyond the MFP 
the jet overcollimates, causing both density and pressure to 
increase again. 

Fig. [3] shows the geometry, the energetics, and the 
Alfvenic Mach number of the jet. The reference field line, 
plotted in the upper left panel, shows that for most of the 
expansion the jet is parabolic with the height of the jet ap- 
proximately the radius to the power of 4/3. At its highest 
point beyond the MFP, the jet has overcollimated to only 
2 per cent of the maximum width, attained just before the 
MFP. The main difference between this solution and solu- 
tion c from PMM10 is that the height-to-width ratio has 
increased by a factor of 2.5, which means it has a narrower 
opening angle comparatively. 

The upper right panel of Fig. [3] shows the partition of 
energy of the jet. After a small drop, the kinetic energy first 
increases at the expense of the thermal energy, signifying 
initial bulk acceleration due to a gas pressure gradient. After 
the flow has cooled, magnetic acceleration takes over and 
continues to accelerate the flow, also after overcollimation. 
Only just before the flow hits the axis does the Lorentz factor 
decrease again with a corresponding increase in the magnetic 
field strength. 

The lower left panel of Fig. [3] shows the opening half- 
angle (ji/2 — tp) of the flow (also compare with Fig.0 and the 
causal connection opening angle. Both angles are measured 
from the axis of symmetry. The opening half-angle shows 
that after a relatively cylindrical start, the flow widens, be- 
fore slowly collimating again. The causal connection opening 
angle is the equivalent of a sonic Mach cone. The flow cannot 
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Figure 1. The reference solution showing both an MFP and 
MSP. The dotted/red and dashed/blue lines show the numer- 
ator and denominator of the wind equation respectively, while 
the solid/black line shows their ratio. The latter determines 
the total bulk acceleration of M 2 with respect to polar angle 
9. The vertical line shows the location of the AP. The verti- 
cal axis is a 'scaled logarithm' of the plotted parameters, i.e., 
y = sign(x') log 10 [l + abs(x)/10 — 12 ] to clearly show the variables 
over many orders of magnitude. At low 8 a small change in angle 
corresponds to a vast change in height, also contributing to the 
near vertical change of sign of the bulk acceleration. (A colour 
version of this figure is available in the online journal.) 



influence anything outside of the forward pointing cone with 
this angle. As the velocity increases, this cone becomes more 
narrow. This plot is almost identical to that for solution c. 

The lower right panel of Fig.[3]shows the Alfvenic Mach 
number, the velocity in units of the Alfven speed. After a 
long initial rise, shortly after overcollimation it decreases 
again. This decrease relates to the leftmost part of Fig. [T] 
where the numerator crosses zero again, without a corre- 
sponding crossing of the denominator, causing a decrease 
of the Mach number. Because the Mach number is defined 
in terms of the Alfven speed, it is not the overall velocity 
that is decreasing, as there is no corresponding decrease in 
the Lorentz factor. It is rather the magnetic field strength 
that increases due to the overcollimation, leading to a higher 
Alfven speed. The plot of the Mach number is very similar 
to that for solution c, but there is a twenty-fold increase of 
M. Note that in PMM10 we plotted the square of the Mach 
number. 
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Figure 2. Various physical quantities that result from integrating the wind equation depicted in Fig. [T] plotted against polar angle 
9. The MFP, AP and MSP are marked. Panel (a) shows the Lorentz factor (7), the poloidal (V p ) and toroidal velocity (V^) of our 
canonical solution. Please note that the toroidal velocity starts out positive, turns negative just beyond the AP (the dotted line) and 
then becomes positive again later in the outflow. Panel (b) shows the electric field energy (E 2 /8tt), the toroidal magnetic field energy 
(B^/87r), the poloidal magnetic field energy (Bp/Sir), the density (pgc 2 ) and the pressure (P). These have all been divided by the scaling 
factor B 2 a F ~ 2 , where Bo is a reference magnetic field strength and a is the square of the cylindrical distance of the AP in terms of a 
reference length (a = ro^/rojj), as defined in VK03. Although the square is plotted, B^ is negative everywhere in this model. (A colour 
version of this figure is available in the online journal.) 




Figure 3. Various physical quantities that result from integrating the wind equation depicted in Fig. [T] now plotted against cylindrical 
radius zu instead of 9. These plots are similar to those in fig. 4 of PMM10. Panel (a) shows the geometry of the field line. Panel (b) 
shows the magnetic energy (S = — zuQB^/'^i j\c 2 ), the thermal energy including the rest mass (£), and the kinetic energy ([7 — 1]£). Panel 
(c) shows the opening half-angle of the outflow (n/2 — ip) and the causal connection opening angle (arcsinp/7]). Panel (d) shows the 
cylindrical radius in units of the 'light cylinder' radius (x) and the Alfvenic Mach number (M). The MSP, AP and MFP (from left to 
right) are indicated for each quantity. (A colour version of this figure is available in the online journal.) 
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Figure 4. Examples of how the heights (in gravitational radii) of the solutions' three singular points change as each of the four principal 
free parameters is changed: solid/red line: MFP; dashed/blue line: AP; dotted/black line: MSP. The initial solution (see Table [3j is 
indicated by the vertical black line. The parameters varied are as follows: panel (a): poloidal spherical angle of the Alfven point (8\); 
panel (b): the angle the field line makes with the disc at the Alfven point (V>a); panel (c): magnetisation (o"m); an< i panel (d): black 
hole mass (AT) in units of Mq. The horizontal axis in each of these plots is linear, not logarithmic. Note the dramatic excursions in the 
height of the MFP (3 orders of magnitude) with only modest changes in the free parameters. Note also that the MSP and AP heights 
are still outside the black hole horizon (see Fig. [5}; i.e., even though the height of the MSP < 2 r g , the MSP spherical radius for these 
solutions remains outside the black hole horizon (r > 2 r g ). (A colour version of this figure is available in the online journal.) 



3.2 An initial exploration of parameter space 

In this section we describe our exploration of parameter 
space. As a starting point we pick our first solution. We 
change the value of a single parameter (8 a, if 1 a, <jm, or Ai) 
until it is no longer possible to obtain a singular solution 
by fitting x\, q, and pa- The effect of changing these four 
parameters on the height of the MFP, AP and MSP is given 
in Fig. [4] 

This plot also shows a physical reason why the solutions 
do not fill up the whole of parameter space. For high values 
of iJja, and for low values of <tm and Ai, the MSP approaches 
the AP. Beyond the location where they coincide, no further 
solution is possible. 

If 8a and tpA are changed together so their sum is 
roughly constant, the range of solutions spans approximately 
45° . If they are changed separately, the range is significantly 
smaller. <tm has the largest spread for this case, ranging from 
0.4 to above 60, monotonically increasing the height of the 
MFP. M has the range 4-12 Mq, spanning most of the 
mass range of astrophysical black holes for an AP distance 
of zd a ~ 74 km. 

Fig. [4] gives the bounds for our solutions. Since we con- 
nect the start of the particle acceleration region with the 
MFP, we are mainly interested in its location within the jet. 
In the solutions found so far, the height of the MFP ranges 
from 10 4 - 10 8 r g . Because this work is by no means an 



exhaustive exploration of parameter space, this range pro- 
vides only a preliminary indication. Nevertheless, it already 
allows us to model a variety of sources. The best way to 
change the height of the MFP seems to be adjusting 8a- The 
height monotonically decreases while increasing this param- 
eter. Another possibility is decreasing <tm- "0a and Ai seem 
to have little effect over a wide range of values. The MSP 
should lie in or near the corona, thus locations too far away 
and too close to the black hole should be avoided. The height 
of the MSP ranges from 2 - 10 r s . This lower value is very 
close to the Schwarzschild radius and justifies the inclusion 
of the pseudo-Newtonian potential. The parameter Ai has 
the largest effect on the location of the MSP. 



The effect of changing Ai on the geometry of the field 
lines is shown in Fig. [5] This plot shows only the single 
reference field line of a solution. It is clear the height of the 
MSP can be smaller than the Schwarzschild radius, but the 
spherical radius can not. This geometry explains why the 
height of the MSP can be smaller than 2 r g in Fig. [4] As 
shown by the lower right panel of this figure, the solutions 
with a higher black hole mass collimate first, and reach a 
lower total height (in r g ). 
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Figure 5. The effect of changing M on the field line geometry. The solid black/line is our canonical solution in Table|3] M ranges from 
7-10 Mq in steps of 0.5 Mq. In panel (a), the MSP and AP are indicated by the black dots. The arrows point to the AP. In panel 
(b) the MFP is indicated by black dots. Note that the height of the MSP (indicated at the lower end of the line) can be smaller than 
the Schwarzschild radius, but its spherical radius always remains outside the horizon. (A colour version of this figure is available in the 
online journal.) 



3.3 Self-similarity 

By including gravity in the relativistic equations of MHD, 
the assumption of self-similarity along the field line is bro- 
ken. Near the black hole, the flow is non-relativistic and the 
wind equation reduces to the form of VTST00, which is self- 
similar. Far away from the centre, the gravitational energy is 
negligible and the wind equation assumes the form of VK03, 
which is also self-similar. This combination means that along 
the integrated field line, the flow transitions from one form 
of self-similarity to another. It should be noted that the ref- 
erence field line we integrate will always be continuous. 

Since one field line does not constitute a jet, we would 
like to be able to describe the flow in a region around our 
reference field line. Changing field lines is done by varying 
the parameter a = zu\/wq, where vja is the cylindrical 
radius of the Alfven point for a specific field line, and ujq is 
a reference length. By choosing this reference length to be 
the cylindrical radius of the Alfven point for our reference 
field line, without loss of generality that field line has a = 1. 

Changing ot to select different field lines then becomes 
equal to changing vja or M. Since this parameter is one of 
our model parameters, the other model parameters have to 
be fit in order to obtain a solutions passing through all three 
singular points. Although the best result would be obtained 
by changing all parameters simultaneously, for simplicity we 
limit ourselves to two. We are free to choose any two param- 
eters, and when using x\ and q the field lines do not cross, 
which is required for self-similarity. Because 6 a and i/ja are 
therefore constant for all field lines in a region around our 
reference field line, at the AP self-similarity is exact. We 
will define this region, where the field lines do not cross and 



satisfy self-similarity to within a specified error, as a flux 
tube. 

To show that the non-crossing of field lines holds for a 
large part of parameter space, we have performed this test 
at different parameter values. Fig. [5] shows the deviations 
from self-similarity of this solution by dividing several field 
lines around a reference field line by this central field line. 
For perfect self-similarity this ratio should be a constant. In 
our case self-similarity is maintained for the majority of the 
jet. Only very close to the MSP do the deviations become 
more pronounced, due to gravity. These deviations are also 
demonstrated by the MSP occurring at slightly different an- 
gles for different values of the black hole mass M. (see Fig. 
O. If we allow a maximum deviation of 10 per cent, which in 
this case will occur at the MSP since we will terminate our 
solutions there, the width of a flux tube would be about 0.30 
of the radius of the central field line, which is a significant 
fraction. 

Beyond the AP, deviations can be caused by changes 
in the fitting parameters to obtain singular solutions for the 
different field lines. The bigger the differences in these pa- 
rameters, the greater the deviations and the narrower the 
resultant flux tube. The required changes in the param- 
eters depend on the location in parameter space, and so 
the allowed width depends on the parameter values chosen. 
For the parameters of Fig. [6] the deviations beyond the AP 
are very small. While our equations do not strictly satisfy 
self-similarity, we have shown that it holds very well and 
that determining the width of the flux tube where it does is 
straightforward. 
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Figure 6. The ratios of the radial size of several field lines with 
regularly increasing vj\ fitted with x 2 ^ and q. The bottom line has 
a = 0.64 and the top line has a = 1.6 (where a = w^/m^). The 
horizontal dashed lines show the values for exact self-similarity. 
The deviations from self-similarity only get pronounced near the 
MSP. The shaded region shows the width of the flux tube for 
which the combined deviation is smaller than 10 per cent beyond 
the MSP. The parameters of the reference solution are x 2 ^ = 0.01, 
o- M = 0.01, q = 0.01, ro A = 5.08 r g , M = 10 M©, 9 A = 0.952915, 
</>A = 89.1119. The Lorcntz factor at the MFP is just above 11, 
showing this solution is rclativistic. (A colour version of this figure 
is available in the online journal.) 



4 DISCUSSION 

By including gravity in the MHD equations, we are able 
to extend earlier solutions that crossed the MFP and AP 
to also go through the MSP, allowing us to describe a jet 
that crosses all three singular points. Because of this lower 
boundary condition, we now have a reliable description of 
the jet below the AP. This description gives a smooth so- 
lution from very near the central object out to the point of 
overcollimation, using a single partial differential equation 
to describe all regimes. This approach produces a workable 
physical model, allowing us to tie the jet properties to the 
conditions at the base, and providing a self-consistent deter- 
mination of the start of the particle acceleration region. 

At the same time, however, the addition of gravity also 
violates the conditions for self-similarity. The reason why a 
single self-similar relativistic flow equation, with gravity, has 
not been derived is because relativistic flow has one scaling 
with radius, while gravity has another. Our C°°-continuous 
bridging method has not changed that situation: while de- 
scribable with a single, continuous equation, our solutions 
in the relativistic part of the flow without gravity will have 
different dependencies on the radius parameter a than in 
the non-relativistic part near the black hole that includes 
gravity. That is, at least one term in our wind equation will 
have a dependency on the radius parameter a. 

In effect our modified wind equation creates two regions 
with self-similar geometry, but with different self-similar de- 
pendencies in each region. Because our focus is on the bulk 
acceleration and collimation of jets in relativistic sources, 
our approach to this issue is to choose the self-similarity 
of the relativistic VK03 equations and restrict the different 
dependency on a to the gravity term only. That is, our solu- 
tions will not be strictly self-similar in the low-speed part of 



the flow with gravity (the VTST00 regime), but they will be 
self-similar far from the black hole in the VK03 regime (see 
Fig- El)- This highlights another advantage of this approach: 
since the gravity term is an algebraic one only (not involving 
any derivatives of the flow parameters with respect to either 
radius or polar angle), all the physical radial and angular 
dependence of the original equations will be preserved. 

The different self-similarities within a single solution is 
not unique to our approach; it will be true of any method 
that attempts to bridge self-similar solutions with and with- 
out gravity. In this regard, interpretation of the solutions will 
be an important aspect of this study. There are two main 
issues to consider. 

First, while accretion discs and jets may have a self- 
similar character, both simulations and observations show 
that their activity is concentrated in specific regions. For 
example, much of the accretion, and therefore jet, power is 
generated near the disc inner edge. And features like the one 
we seek (a strong collimation shock in the flow) will occur 
primarily at a specific point in the jet, and be observed as 
such, rather than being spread over a large range in radius. 
Therefore, we shall concentrate on only one bundle ('flux 
tube') of magnetic field and stream lines and assume that 
this flux tube is anchored near the disc inner edge and that 
it passes through the most important features of the jet. 
This will limit the range in a in which we are interested. 

Second, in order that the solutions remain reasonably 
self-similar-like, we shall choose solutions in which the field 
lines do not cross within that flux tube. By choosing a spe- 
cific combination of the free parameters to fit a solution 
with three singular points, the field lines around our refer- 
ence field line behave well. We have found that selecting x\ 
and q as fitting parameters for the singular points produces 
results that satisfy self-similarity well within a flux tube of 
finite width. 

With our method, therefore, we will restrict our solu- 
tions to a limited, and conservative, region of parameter 
space. The flow in our solutions must remain reasonably 
non-relativistic in the VTST00 part of the flow (when grav- 
ity is important), and we must consider only flux tubes of 
narrow enough width that field lines satisfy self-similarity 
to within a specified error in this region with gravity. 

Fig. [7] can be used to see whether the assumptions of 
the two regimes are not violated for a particular solution. 
The kinetic energy (7 — 1) should be negligible near the 
black hole, while the gravitational potential energy (%^-) 
should be significant, satisfying the conditions of the non- 
relativistic regime and demonstrating the importance of 
gravity. Around the Alfven point the kinetic energy over- 
takes, without becoming relativistic. This is the bridging re- 
gion where gravity becomes unimportant. Only beyond the 
Alfven point does the flow become truly relativistic, while 
gravity becomes negligible, satisfying the conditions of the 
relativistic regime. 

In this work we have only made a cursory exploration of 
the full parameter space, and it seems likely that part of that 
space will allow for solutions with smaller MFP values. The 
location of the MSP falls within the range of a radiatively- 
inefficient accretion flow (RIAF) disc model. Also, the ve- 
locity at the MSP, around c/3, is very reasonable. These 
features make it plausible that we will be able to match the 
conditions to an accretion flow as well as to the size of the 
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jet b ase as indicated by fits to the spectrum (|Markoff et al.l 
120051 ). allowing us to determine the location of the shock 
region from the conditions at the base of the jet. 

It is important to note that since our assumption of non- 
relativistic flow (when gravity is important) eliminates the 
part of parameter space that produces relativistic flows close 
to the black hole, we exclude solutions that may relate to 
physical sources. Future work would include a more general 
model that allows relativistic flow in the gravitational field. 
However, the model in this paper provides a good description 
of the jet structure near the MSP when the flow is non- 
relativistic. 

We should also note that the gravitational field we chose 
to include here is, at best, a pseudo-Newtonian (or pseudo- 
Schwarzschild) one. Our solutions, therefore, do not take 
into account any metric rotation that would occur near a 
Kerr black hole, for example. However, our method might 
be able to hint at the presence of a Kerr hole, if, for example, 
after fitting to broad-band data, we find that the footpoint 
of the field line may be significantly less than the ISCO of 
a Schwarzschild black hole (<g 6r g ) or that the required 
magnetisation is significantly greater than would be typical 
at the ISCO of a normal Keplerian disc. 

Another issue neglected by our computations (and, in- 
deed, by all steady state analyses) is the effect of insta- 
bilities on the accelerating jet flow. We will briefly discuss 
fluid and MHD instabilities in the weak-field limit and then 
MHD instabilities in the strong-field limit as well. Figure 
2b compares the relative strengths of fluid dynamical forces 
(i.e., pressure P) with electromagnetic forces (e.g., (B 2 + 
E 2 )/8n). We see that the former dominate only below the 
Alfven point, in the slow magnetosonic region, so it is only 
there where we expect Kelvin-Helmholtz (KHI) or magneto- 
rotational instabilities (MRI) possibly to be important. In- 
deed, the MSP itself may lie in the atmosphere of a turbu- 
lent accretion disk, giving rise to a somewhat more complex 
structure in the sub-Alfvenic region than that assumed here. 
Nevertheless, recent two- and three-dimensional simulations 
of relativistic jet s launched from turbulent acc retion disks 
ljMcKinnevll2006l : iMcKinnev fc Blandfordl l2009l l show that 
MHD jets similar to those investigated here (including a 
well-defined classical slow point) do indeed form and are 
not disrupted by weak-field instabilities near the accretion 
disk. Above this point, like other outflowing magnetospheres 
(solar, pulsar, etc.), the dynamical forces in the flow are so 
dominated by electromagnetic forces (eventually by many 
orders of magnitude) that any weak-field instabilities will 
be strongly suppressed, or at least unimportant in affecting 
the kinematics of the accelerating jet. 

However, strong-field instabilities (in particular, the 
current- driven helical kink [CDI]) n eed to be looked at more 
closely. iNakamura fc Meierl l)2004h . for example, showed 
that in non-relativistic jet flow rotation velocities well 
above the Alfven speed can suppress the helical kink, 
as can a steep external pressure gradient. They specu- 
lated that relativistic flow would further reduce the de- 
velopment of kinks, but a similar detailed relativistic 
study has not yet been performed. Again, as one of 
the few relativistic three-dim ensional MHD simulations, 
IMcKinnev fc Blandfordl (|2009l ) can give some insight into 
the behavior of real MHD jets in the strong-field region. 
While the accelerating jet appears to be affected somewhat 
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Figure 7. Comparison of the gravitational potential energy and 
the kinetic energy along the jet of the bottom solution detailed 
in Table [3] With this plot we can select solutions where the grav- 
itational energy QA4 / (c 2 r) is important close to the black hole, 
whereas the kinetic energy (7 — 1) only becomes relativistic fur- 
ther along the jet. (A colour version of this figure is available in 
the online journal.) 

by helical kinking in that study, it does not appear to be 
disrupted by the CDI. It remains a viable jet well beyond 
the classical fast point (where the rotational speed should be 
significantly super- Alfvenic. However, no three-dimensional 
simulation so far has followed the acceleration out to the 
MFP region, let alone for long model times to achieve a 
quasi-steady state. And, detailed numerical parameter stud- 
ies are even further in the future. So, the question of whether 
MHD jet acceleration can be largely stable to strong-field 
MHD instabilities is still an open question. 



5 CONCLUSIONS 

We have shown that it is possible to extend the time inde- 
pendent, semi-analytic solutions of PMM10 to solutions that 
include gravity using a C°°-continuous bridging method, 
connecting two valid regimes with a smooth transition re- 
gion. We also extend the MHD wind /jet equation to include 
a pseudo-Newtonian potential for gravity. For the first time 
these solutions cross all three singular points and describe 
a relativistic jet from the vicinity of a central black hole to 
the point where the flow hits the axis of symmetry. Although 
these solutions do not satisfy strict self-similarity, the devia- 
tions within a flux tube of a certain width can be estimated 
and controlled, and thus a physically relevant model can be 
constructed. 

We have discovered a solution crossing all three singular 
points with parameters very close to one from our previous 
paper. This new solution shows the relative ease with which 
solutions can be found and allows for a comparison. The 
main differences are the higher distance of the MFP and a 
significantly higher velocity, both owing to a higher bulk ac- 
celeration at the AP. Otherwise the two solutions look very 
similar, showing that gravity has little influence beyond the 
region close to the central object. By allowing the creation 
of an MSP, including gravity enables us to have a reliable 
description of the flow below the AP. 

We have made a cursory examination of parameter 
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space by changing one parameter of a single solution at a 
time. This exploration gives a preliminary indication of the 
extent of the solution space, while at the same time show- 
ing the effects it has on the solutions. In particular we have 
found that the poloidal spherical angle of the AP has the 
most significant effect on the height of the MFP. On the 
other hand, the ratio of the AP to the light cylinder radius 
has the greatest effect on the height of the MSP. 

The multi-dimensional solution space is constrained by 
physical and mathematical conditions. Changing only one 
parameter at a time, when the light cylinder radius ap- 
proaches the AP, the AP moves outwards, the temperature 
is increased, or the magnetisation decreased, the jets become 
wider and shorter, with the MSP approaching the AP. By in- 
creasing the cylindrical radius of the AP, the effects of grav- 
ity decrease. Eventually this decrease leads to the regime 
where an MSP cannot be created anymore, approaching the 
situation without gravity. Despite all these constraints, the 
solution space allows a wide range of parameter values to 
be chosen, translating into a wide range of properties, like 
the location of the MSP and MFP, the velocity, magnetic 
energy, density and pressure of the flow. 

By matching conditions at the MSP to an accretion flow 
model, the mathematical parameters considered free in this 
work will be tied to the conditions at the base of the jet and 
subsumed into the model, with the wide range of properties 
in the solution space assuring a good fit. Consequently, the 
height of the MFP, and the start of the particle acceleration 
region, will be uniquely determined by the conditions close 
to the central object, providing a self-consistent connection. 
After integration into a model that can determine the spec- 
tral emission of a jet solution, we will be able to ascertain 
the conditions that best describe the overall spectrum of a 
given black hole system. In future work we will explore how 
well this model succeeds in predicting the correct location of 
the optically thick-to-thin break observed in the broadband 
spectra of compact jets, hopefully with new insights about 
the physics of jet launching conditions. 
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Substituting these into equation (17) of VTSTOO pro- 
duces: 



GM 

m*V 2 



GM P?x\ (1 - M 2 - x\) 2 
z 2 vj A F 2 a 2 1 (1 - M 2 - x 2 ) 2 ' 



The full gravity term then becomes: 

_ n 2 sin(fl) _ GM n 2 xj (1 — M 2 — x\) 2 sin(fl) 
G c 2 F 2 a 2 A (1 — M 2 — x 2 ) 2 vj k G ' 



(Al) 



(A2) 



In this expression (wxG)j sin(#) is equal to the spherical 
radius, r. To include a pseudo-Newtonian potential, we re- 
place 1/r by l/(r — rs), where rs is the Schwarzschild radius 
(= 2GM/c 2 ) and simplify: 



k 2 sin(0) _ n 2 x\ (1 — M 2 — x 2 A ) 2 
G ~~ ~^jT (1 — M 2 — x 2 ) 2 



GM sin(6>) 



(A3) 



APPENDIX B: EQUATIONS FOR THE INITIAL 
PARAMETER VALUES 

The Alfven regularity condition is given by evaluating the 
wind equation at the AP: 

FVj^l-ajKaA + l^sin^A) 

/i 2 COS 2 (# A + IpA) 

2 r^(F-2)(g A ^i )( i-,i) s . n(gA) 



+ 2 cos(V>a) sin(0 A + iPa) 



CTA + l 
0"A 

+ *EV*1 [(F -i )( i -. a) -i] 



= [x A - cta(1 - x A )] - (F - l)cr A (l - x A 

- 2 



r 

2 

x A 



1 



(Bl) 



6 

C 2 Z&A 



GM sin(0 A ) 



- 2 



We obtain fi by evaluating the Bernoulli equation at the 
AP: 



(cta + iy 



- [x\ - GA (1 " x\)] 2 

F 2 <Jm (l-z A ) 2 sin 2 (^ 



+ 



2 fl 



x\ cos 2 (6 i a + Va) 



(B2) 



APPENDIX A: DERIVATION OF THE 
GRAVITY TERM 

To obtain a relativistic form of gravity, we need the following 
conversions: zn t in VTSTOO is equal to vjq in VK03 and V* 



similarly corresponds to — 



'*Fcr M c 

-rf* A 



APPENDIX C: GRAVITY IN THE BERNOULLI 
EQUATION 

The only terms that include the effects of gravity in 
the Bernoulli equation given by (A3) in VTSTOO are the 
Bernoulli constant and the gravity term. Dividing out a com- 
mon factor of 2, we have: 



£ + 



k 2 sin(fl) 
G ' 



(CI) 
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At the AP G = 1 and 6 = 9a, reducing these terms to: 

e + K 2 sin(6lA). (C2) 

If we write out e using equation (19) of VTSTOO, G = 
zu/zua, and equation (2.7a) of BP82, we obtain (please note 
that to keep to the notation of the preceding two papers, 
here the subscript means the value at the base of the 
outflow, not the reference values as used in VK03. Unfor- 
tunately, ro as used in BP82, is vjq as used in VTSTOO): 



Go QM /wo vjq/kja 



GM/uja ' 



(C3) 



Since e, the specific energy, is a constant of motion, we can 
evaluate it at the AP, just like the Bernoulli equation in 
VK03 (equation (2.2) of BP82 with the gravitational poten- 
tial expanded): 



Vl , , QM 



e = e A = -y + fiA - 



Quja B, 



0,A 



/'A 



*A,A 



(C4) 



Here ta is the spherical radius of the AP. Substituting (|C4[) 
and (IC3l) into (IU21). we obtain: 



VJA 

QM 



ZUA , • la \ 
hsm(6» A ) 



*a,a J r A 

(C5) 

However, zoa/ta — sin(#), so the last two terms cancel. The 
k term cancels the QM factor (equation IA1|I , leaving no 
terms that depend on gravity. Thus there is no dependence 
on gravity in the Bernoulli equation in the non-relativistic 
regime with gravity at the AP and consequently there is no 
gravity term in the relativistic Bernoulli equation at the AP 
(equation IB2[) . 
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